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Abstract 

Two-qubit X-matrices have been the subject of considerable recent attention, as they lend them¬ 
selves more readily to analytical investigations than two-qubit density matrices of arbitrary na¬ 
ture. Here, we maximally exploit this relative ease of analysis to formally derive an exhaustive 
collection of results pertaining to the separability probabilities of generalized two-qubit X-matrices 
endowed with Hilbert-Schmidt and, more broadly, induced measures. Further, the analytical re¬ 
sults obtained exhibit interesting parallels to corresponding earlier (but, contrastingly, not yet fully 
rigorous) results for general 2-qubit states-deduced on the basis of determinantal moment formu¬ 
las. Geometric interpretations can be given to arbitrary positive values of the random-matrix 
Dyson-index-like parameter a employed. 
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I. INTRODUCTION 


In previous work l|, 2], the authors investigated the separability probabilities of 4 x 4 


density matrices (p) with density approximation techniques [3j based on the moments of the 
determinant (|p PT |) of the partial transpose-the nonnegativity of which is necessary and 
sufficient for separability jd, 5]. The moment formulas employed involve a Dyson-index-like 
parameter a which specializes to real, complex, quaternionic densities for a — |, 1 and 2, 
respectively. As yet, the validity of the moment formulas is only a conjecture, but there is 
strong evidence that they are correct. By use of extraordinarily large (5 x 10 11 trials) Monte 
Carlo simulations, Fei and Joynt [6] obtained separability probability values agreeing to 4 


• Blab, 


64 , 33 323 , respectively. More__generally 




decimal places with our results (cf. Id, |8|), that is It, and 26 

sti "’ we have al8 ° been mvest,gatmg the general separabmty probabllity que8tl011 fl when 

th powers of the determinant (|p|) times the Hilbert-Schmidt (HS) measure [10j, [llj-that is, 


induced measure 


3,Q, 


13]-are used as the probability measure on the convex set of density 


matrices. Due to the extreme intractability of performing direct probability calculations by 
use of high-dimensional integration, these studies continue to be based on moment formulas 
and sophisticated numerical analysis (density approximation) techniques B to determine 
probabilities. 

Recently, interesting results have been obtained for the subset of so-called X-matrices, 
which are a form of toy model for density matrices 7, 14,il5|. The present study exploits their 
simpler structure to directly, and now fully rigorously, compute the separability probabilities 


for the family of induced measures (including the particular [k = 0] Hilber 


Although the results obtained are numerically quite different from those 16] holding for the 


-Schmidt case). 


full, unrestricted matrices, they do exhibit a strong qualitative similarity. This can be taken 


as evidence that our earlier work 




16] has been well directed. We note that in the 


X-matrix case, it has been possible to give geometric interpretations to arbitrary positive 
values of the Dyson-index-like parameter a , not just integer and half-integer values, such as 
1 and 2 as described above (cf. 17]), and perhaps this idea can be extended to unrestricted 
sets of density matrices. 

We start with the basic definitions. A 4 x 4 X-density matrix £ is positive semi-definite 
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with Tr (£) = 1 and has the form 


61 0 o fu 

£ _ 0 £22 £23 0 

0 £23 £33 0 

<^14 0 0 £44 _ 

the defining conditions are equivalent to > 0 for 1 < j < 4, £ 14 , £23 € C, Y2j=i£jj = 1; 
|^i 4 1 2 < ^ 11^44 and |^ 23 1 2 < ^ 22 ^ 33 - The partial transpose of £ is denoted by £ PT and 


£ 


PT 


£11 0 0 £ 23 

0 £22 £14 0 

0 61 £33 0 

£23 0 0 £44 


Thus det£ pr = (£n £44 — |£ 23 1 2 ) (£ 22£33 — |£i 4 | 2 ) and the density matrix £ is separable if and 
only if det £ pp > 0. 

We induce a measure on the set X of X-density matrices from the measure 

4 

JJ d i j3 (r 5 r 6 ) 2 “ _1 dr 5 dr 6 d6> 5 dl9 6 , 
o=l 


on the cone of all positive semi-definite X-matrices, where £14 = r 5 e 10s ,£23 = r 6 e l6,fi (r 3 > 
0, — 7 t < 6j < 7 r, j = 5,6). The details of normalizing the measure to be a probability 
measure appear later. We will show for any a > 0 that 


m 


Pr {det £ pr > 0} = p (a) := 
particular, p 


T (a + |) T (a + §) 


p(1) - I ^ 7 


^r^ + DPta + ljr^ + f)’ 
2 


The value p (1) = | is also found in 


Milz and Strunz [7|, eq. (22)], where they employed a quite different, interesting analytical 
framework, in which the principal variable of interest is the radial location in the Bloch ball 
of the reduced density matrix. 

Further, we find Pr {det£ pp > 0} when the measure is multiplied by (det£) fc for some 
fixed k > 0 and a = 1,2,3,.... For example when a — 1, Pr {det £ PT > 0} = 
1 — F ( 2 + 2 )r{- 3 4 A-+ 6 ) • Tor a = 2, 3,.. .we find that Pr {det£ PT < 0} is a product of a ratio 
of gamma functions with a polynomial in k with positive coefficients of deg ree 2 a — 3. This 


combination is similar to our results for the full 4x4 density matrices 


161 ]- 
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Section 2 contains the construction of the coordinate system used to set up the relevant 
def ini te integrals as well as the expression of the normalized measure in this system, and 
explains the relation of the values a = |, 1 , 2 to real, complex, and quaternionic matrices. 

In Section 3 the computation of Pr {det£ PT > 0} for any a > 0 is carried out. The 
computation starts with a five-fold iterated integral which is reduced to an integral of hy¬ 
pergeometric type and finally to a classical hypergeometric summation formula. 

The measure of (det^ type is studied in Section 4. The calculation of Pr |det£ PT < 0} 
is carried out only for integer values of a for reasons of technical difficulty. As will be 
seen, even this integer case is quite complicated. Necessary integral formulas are derived in 
Section |V] 


II. THE COORDINATE SYSTEM AND THE MEASURES 


We construct a family of probability measures on A, with a parameter a, which 
agree with the normalized measures induced by the Hilbert-Schmidt measures on 
M 4 (M), M 4 (C), M 4 (H) for a = A, 1 ,2 respectively (H denotes the quaternions). As in 
our previous studies 


ilii 

EL 


16 


18], we use the Cholesky decomposition to define the measures. 


For x G M > 0 x C 2 let 


xi 0 0 X5 

0 x 2 %6 0 

0 0 £3 0 

0 0 0 x 4 


f = C*C = 


x\ 0 0 aqxs 

0 x| x 2 x 6 0 

0 x 2 xci xl + |x 6 | 2 0 

X 1 X 5 0 0 x\ + |x 5 | 2 


eventually we will impose the restriction x j + l x 5| 2 + 1^61 2 = 1 so that trC = 1, but 
first we work with arbitrary positive-definite matrices. 

To motivate the definition of the measures, we describe a simple model for the 2x2 
block. Consider the map defined on the subset H : = x M m of M 2+m (for some fixed 
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m = 1 , 2 ,...) 

</>: (x 1 ,x 2 ,yi, ■■■ ,y m ) H- (x 2 ,^ + \y\ 2 ,x 1 y 1 ,... ,aq?/ m ) ; 

The Euclidean measure on H can be expressed as t\ 1 ^ 2 t. 2 ] ' 2 t^ 2 1 dtidt 2 dt 3 duj {y') where 
t\ = x ' 2 , t ‘2 = x\ and t 3 = \y\ 2 ,y = \y\y' and du {y') is the surface measure on the unit 
sphere { y ' G R m : \y'\ 2 = l}. The point ( t 4 ,t 2 ,t 3 ) G R+. The Jacobian of 0 equals 4xJ + ”\t 2 . 
The image of the measure on H under 0 is (constants are discarded) 

t( n / 2 f™/ 2 1 dt\dt 2 dt 3 du (y f ). 


Now adjoin another copy of H and the map (a direct sum) and relabel to arrive at 
0: (xi,x 2 ,x 3 ,x 4 , 2 / (5) , 2 / (6) ) G R+ x R 2m (xl, xj, xj + \y (6) f , x\ + |y (5) | 2 , x x y^\ x 2 y {6) 


and the measure 


t™l 2 t™/ 2 t™/ 2 1 t 1 ^ 2 1 dtidt 2 dt. 3 dt 4 dt 5 dt 6 du (y^') du (y^') . 

In the cases m = 1,2,4 this construction can be interpreted in terms of the Cholesky 
decomposition of a 4 x 4 positive-definite X-matrix over R, C, HI respectively. In this situation 
det£ Pi = (t 4 t 4 + t\t 5 — t 2 t e ) (t 2 t 3 + t 2 t 6 — t\t 5 ) and the du factor does not enter into the 
calculation, and so is replaced by 1. The same result is obtained if £i 4 , £23 are replaced 
by (fit 5 ) 1 / 2 e ie5 , (t 2 t 6 ) 1 / 2 e ie6 respectively (with — n < 0 5 ,9 6 < n) and du (y^ 1 ) du (y^') is 
replaced by (^r ) 2 d0 5 d0 6 . To sum up this discussion, the generic 4x4 positive-definite 
complex X-matrix is 

t 4 0 0 (tit 5 ) 1/2 e i05 

0 t- 2 (t 2 te) 1/2 e 106 0 

0 (t 2 t 6 ) 1/2 e~ 106 t 3 + t 6 0 

(fits) ^ e i05 0 0 t 4 -\- 1 3 

and the parametrized measure is 

2 f 5 _ 1 t$~ 1 dt 1 dt 2 df 3 df 4 df 5 df 6 d< 9 5 d 0 6 ; 


which has geometric interpretations when a — y, with m = 1, 2,4. The last step is to induce 
this measure on {£ : Tr£ = 1}, that is, on the unit simplex T 5 in R 5 (t 2 — 1 —1 4 — X/y= 3 tj) 
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and drop dt 2 from the measure. Also since we are concerned with only det £ PT and det £ we 
also drop dO^dde. 

We have arrived at the measure on T 5 

d n a = t*t%t%- 1 t%- l dt 1 dt 3 dt 4 dt 5 dt 6 , 

with normalization constant (a Dirichlet integral) 

T (4a + 4) 
r(a + l) 2 r(a) 2 ’ 

1 . 

With these coordinates 



det £ = t l t 2 t 3 t 4 , 

det £ P1 = (t 4 t 4 + t\t§ — t 2 te) (t 2 t 3 + t 2 tft — tits). 

We introduce the desired coordinate system in two steps. The first step is: 

£n — t\ = - (1 — Si + s 2 ), £ 22 = t 2 = - (si + S3) , 

1 1 

£33 = £3 + h = - (si — s 3 ), £44 = t 4 + t 5 — - (1 — Si — s 2 ), 

1 1 

t$ = -S 4 (1 — Si — S 2 ) , t 6 = -S5 (si — S 3 ) , 

^3 = g (1 — S5) (Si — S 3 ) , t 4 = - (1 — S 4 ) (1 — Si — S 2 ) , 

where 0 < si, s 4 , S5 < 1 and |s 2 | < 1 — si, |S3 1 < Si. The Jacobian (omitting t 2 from the list 
(t) ) is 

= 1 (S ‘_ S3) a _ Sl _ S2) , 

and the measure and det £ PT transform to 


d ^ = 2 4a 2 ((1 - si) 2 - s 2 2 ) a (si - sl) a S4 ^5 1 dsids 2 ds 3 ds 4 ds 5 , 

det i PT = ~j ~0 (((! - s lf - si) - S 5 (s? - -8) ((s 2 - 4) - s 4 ((1 - Si) 2 - si)) . 

Since det £ PT is even in s 2 and s 3 we can restrict to 0 < s 2 < 1 — Si, 0 < s 3 < si (a “quarter” 
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of X, denoted by Xq) and multiply the measure by 4. The second step is to change variables 


s 2 = yj (1— si) 2 '— 4hi, s 3 = \Js\- 4h 2 , (1) 

d(s 2 ,s 3 ) = _4_ 

a(<Ma) V / (l-, 1 ) 2 - 4 h i yi^ 4 ^’ 

det i PJ = (hi - s 5 h 2 ) (62 - s 4 hi), 
det £ = hih 2 (1 - s 4 ) (1 - s 5 ), 

det £ pt - det £ = (hi - h 2 ) (s 5 h 2 - s 4 hi). (2) 


The measure transforms to 


a— 1 „a—l 


dv„ = 


Xa£a a-i a 
°1°2 6 4 6 5 


|(1 - si) 2 - h M \s\ - h 


:dsidhidh 2 ds4ds5, 


0 < Si,S 4 ,S 5 < 1 , 0 < hi < 


1 - Si 


0 < <s 2 < (4 


There is a beta-integral which we will use later, and also to verify the normalization (recall 

B (°> h ) : = ixSSr = fo pa_1 - p ) 6_1 dp ) 


u 


1 0 y/c — U 


du = d 


8+1/2 B f 


V 


« + !,- 


/ o—2—4 a / I \ 2 »1 

dv a = ^ B ( a + 1 , -) J (1 - Si ) 2 " +1 s 2Q+1 dsi (3) 

_ 2~ 2 ~ 4a T (a + l ) 2 7 tT (2a + 2) 2 _ 1_ 

« 2 T (a + |) 2 T (4ct + 4) c a ’ 

by use of the T-duplication formula T (2 u) = J=2 2u-1 T (u) T (u + |). We will also use the 
Pochhammer symbol, (t) 0 = 1, (t) n+1 = (t) n (t + n) for t € C. 

We finish this section by describing the extreme values of det £, det £ PT and det £ PT —det £. 
The maximum value of both det£ and det£ Pi is -W for £ = 1/ (identity matrix). The 
minimum value of both det£ PT and det£ pp — det£ is — for si = S 5 = 0 , S 4 = l,hi = 
4, h 2 = 0 (and other matrices). To maximize det£ PT — det £ = (hi — h 2 ) (ssh 2 — S 4 hi) set 
s 5 = 1 , s 4 = 0 , hi = ) h 2 = (y ) 2 to obtain ys 2 (1 — 2 si) with maximum value ^ a t 

Si = |. These extreme X-state values are identical to those for the full matrices. 
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III. THE COMPUTATION OF Pr{det£ PT > 0} 


The desired probability is the i^-measure of the set 
{(si, Ai, <5 2 , s 4 , s 5 ) : (Ai — S 5 A 2 ) (<5 2 — s 4 Ai) > 0 }, in other words, the definite integral of 
du a over this set. We start the iterated integral with the variables S 4 , s 5 . There are 
apparently two possibilities for (<5i — s 5 <5 2 ) {82 — s^i) > 0 : 


^1 << o _ ^ o. (§1 §2 \ - j _ 


<5i — s 562 < 0 and 82 — S 4 (fi < 0; this implies ^ < S 5 and y < S 4 , but max , s , ^ 
and one of the inequalities contradicts s 4 , S 5 < 1 , except for the trivial case <5i = 
8 2 , s 3 = 1 = s 4 , included in the following case (the products of the two pairs of 
eigenvalues of £ PT are (<5i — S 5 < 5 2 ) and (82 — 54 ^ 1 ), so this demonstrates that £ PT can 
have at most one negative eigenvalue); 

— s 5 8 2 > 0 and 8 2 — s 4 8 1 > 0; equivalent to ^ > s 5 and > s 4 . Thus 8 2 < c>i 
imposes the bounds 0 < S 4 < < 1 and 0 < s 5 < 1 < Ju Similarly 5] < 8 2 

imposes the bounds 0 < S 4 < 1 and 0 < S 5 < < 1. In both cases the integral of 

S 4 _ 1 S 5 ^ 1 ds 4 ds 5 over this region equals 


a i 


( min (|- 


82 


As a side observation, the computation of Pr |det£ PT > det£} = 
Pr {(<5>i — 8 2 ) (s 5^2 — S 4 < 5 i) > 0} starts with integrating S 4 _ 1 S 5 _ 1 ds 4 dss over the region 
8 \ > 5 2 ,0 < s 4 < yS 5 , or the region S 2 > <5i, 0 < s 5 < ^s 4 ] the result in both cases is 


1 

2a 2 


<5 2 <5i 
<5i ’ <52 


mm ( ) | . Since the rest of the probability calculation is the same for both, 


we see that Pr {det£ PT > det£} = | Pr {de t £ PT > 0}, analogously to the general 4x4 

181 - 


situation, as we discussed in earlier work 

By symmetry, it suffices to integrate over {0 < si < |} and double the result. Denote 
a := (^" 2 £1 ) 2 and b := (y) 2 , thus 0 < b < a < |. In the iterated triple integral first integrate 
with respect to <fi. Let 


2 f 1/2 

I a := —z / dsi 


8 a 

°2 


« 2 JO 


- 1/2 


a- 


dsi 


*0 yb — 82 
f b d 8 ‘ 


=d5 2 


'0 


f 62 ( 8 j_ 

0 V^2 

<5 2 § 2a 


V^ dS ' + L (I) 


dhi 


\Jb — 82 0 V 0 , — hi 


dhi + 8 2a 


d<5i 


2 


yja - hi 
















In the {•} expression, the second integral equals 2 yja — S 2 ; for the hrst one, interchange the 
order of integration (over 0 < hi < 62 < b ) to obtain 


rb r2o 


K 


zdtfi 


dh 2 


\J a - hi J 5l y/b - 62 


= 2 ^ 


2a 


6-fr 

a — hi 


Thus, 


I = — 

la .2 


' 1/2 


a z 


dsi / v 


2 a 


b — v / a — v 

a — v V b — v 


dhi. 


dt>. 


This integral can be directly evaluated for any a with 2a G Z + with elementary methods. 
However, this integral is one of a parametrized family which all have closed forms for their 
values. The integrals are denoted by 


I(m,n):= / ds v m < (a — v) 


n ,a-v , ,, . n lb — v 


b — v 


+ {b-v) r 


a — v 


dv. 


The formulas for (m, n ) = (2a, 0) for a > 0 and for m, n — 0,1, 2, 3,... are derived in Section 

El 


From Proposition IV.21 it follows that 
la = —zl (2a, 0) = 


T(2a)T(2a) 


a- ' ■ 2 2+8q T (2a + |) T (2a + |) 

When 2a G Z + the reciprocal of I a is an integer. In fact, 

IT 1 = 


1 — 2(4a + 3) (^°2 


Finally, we obtain 

Pr {det £ PT > 0} = c a I a = 


7T 


T (4a + 4) r(2a)r(2a) 


2 2+8q r (a + l) 2 r (a) 2 r (2a + §)T (2a + |) 
1 r(a + i) 2 r(a + |) 


( 4 ) 


\/2tc r(a + f)r(a + l)r(a+f) 


=:p(a) 


by repeated use of the T-duplication formula. If a = m is an integer, then c a = 

(m - l) 2 (4m + 3)! , . . . , T . , , 4 (m — ~) 2 (4m + 2)! 

-—- and p (a) is rational. If a + ^ — m, then c a = — - — j\ - 

mr - tt 2 ± 

\2 / m 

and p (a) n 2 is rational. The asymptotic formula —-- ~ u a ~ b as u —> oo shows that 

T (u + b) 

p (a) ~ ^ as a —> oo. For example, p (10) = 0.12683... and ^== = 0.12616_ 


V 2 


na 
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A. Minimally degenerate matrices 

In our previous work jji|, we considered the separability probability of a minimally 
degenerate density matrix, and found numerical and formulaic evidence that this probability 
is exactly one half of the unrestricted one. We can prove this relation in the present X-matrix 
case. From det£ = hih 2 (1 — s 4 ) (1 — S5), we see that a necessary condition for minimal 
degeneracy is that one of hi, h 2 , (1 — s 4 ), (1 — s 5 ) vanishes. The possibility hi = 0 or h 2 = 0 
is not minimal, for suppose hi = 0 then (see (P)) s 2 = 1 — «i which implies f 4 = 0 = t 5 , 
reducing the dimension by two. This leaves the cases s 4 = 1 or s 5 = 1 (for example, s 5 = 1 
is equivalent to £ 22 £ 33 — (£ 231 2 = 0) which remove only one dimension (degree of freedom). 
So the minimally degenerate subset of X consists of two almost (i.e. the intersection has 
measure zero) disjoint subsets, one with s 4 = 1 and one with s 5 = 1. It suffices to compute 
the probability for one of these, say S 5 = 1. First we compute the normalizing constant: the 
calculation is the same as in (J3J) except for the factor ^ from J 1 s^dss, thus the constant 
is We see det £ P1 = (hi — h 2 ) (h 2 — s 4 hi) > 0 if hi > h 2 and 0 < s 4 < J^, or hi < h 2 and 
s 4 > ^ > 1; as before the second case is impossible. Thus the first step of integration (for 
s 4 ) yields )- for hi > h 2 and zero otherwise. In contrast to the earlier computation, 

there is no Si -H- 1 — si symmetry so 0 < Si < | and \ < Si < 1 are handled separately. 

If 0 < Si < | (so that b = (y) 2 < = o), then the remaining triple integral is 

W 2 r b a? r° 

A. 


1 

a 


2 

a 


dsi 

fL/2 


0 Vb ~ ^2 

r b 

dsi 


=dh 2 


is 2 


h“ 


\J a - hi 


dhi 


r2 a 


'0 


b- 5. 


-dh 2 . 


If | < Si < 1 (so that b > a > Si > 62 ), then the remaining triple integral is 


1 

a 


dsi 


So 


>1/2 

2 -i 

a 


dsi 


1/2 


Vb=S 2 


r2Q! 

^2 


=dh- 


2 h 1 - 


S? 


\/a — hi 


dhi 


a 


b-S . 


-dh 2 . 


In the latter expression change variables s = 1 — si which interchanges a and b with the 
result 


2 

a 


" 1/2 


ds 


r2a 

^2 


b — So 


dh 2 
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and adding the two parts leads to 


2 

a 





2 

a 


I ( 2 a, 0 ). 


Thus 

Pr {det £ pt > 0} = ——I (2a, 0) = -p (ct), 
a a 2 

see equation ([4]). 


IV. THE MEASURE (det£) fc du a 


Here we compute Pr jdet£ pp < 0} when X is furnished with the normalization of the 
measure (det £) fc dzy*, for a,k e Z + . It appears possible to carry out the calculations for 
fixed integers k and arbitrary a > 0, but with the goal of allowing k as a free parameter, 
technical factors impel us to restrict to integer a. Recall det£ = 8182 (1 — s 4 ) (1 — S 5 ). The 
normalization constant is 


Ca,k 


r (4ct + 4k + 4) 


T (k + a + if T (k + l ) 2 T (a) 2 ' 


and the measure is 


§ a + k S a+ k s a-l s a- 1 ^ ^ ^ 

t(1-si) 2 -5i\As?-5 2 


dsid(fid<f2ds 4 ds5. 


Some experimenting suggests that it is more tractable to compute Pr {det£ PT < 0}. Then, 
the first step is to compute (with <5 0 = min (fj, )) 


a -1 (-1 \k A f 1 a-l ft X A r (a) T (k + 1) f 1 x ( ^ k 

u (1 - u) du / v (1 - v) dv = —— v (1 - v) dv, 

Js 0 T(a + k + l) J 5o 

where {u, n} = {s 4 , s 5 } (depending on whether 81 > 82 )- The second integral is an incomplete 
Beta integral. As discussed above, the restriction a G Z + leads to a feasible calculation. We 
use this simple antiderivative formula: 




a—1 


(\ _ „W+J+1 _ M ,.\k 


3 =0 

v“- 1 (1 - vf 


(k + l ) j+1 

a—1 

dp = 

j =0 


(i-^r 

(1 - a) J 

(k+ 1) 


= v a - 1 (1 - vf , 


- 8 ^ (1 - 8 0 ) k+j+1 . 


3 +1 


(5) 
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From this point on, we work with each term in the sum separately, in the triple integrals 
that remain to be evaluated. At the end, the parts will be summed over j to get the desired 
probability. 

If 5 1 < S 2 then (1 - -<5o ) k+j+1 = 5l a+k ~ l ~ j ( 6 2 - ^) fc+J+1 , otherwise 

6 a + k S a+ k 6 a-l-j ^ _ ^fc+j+1 = ^a+fc-l-j ^ _ § ^k+j+ 1. Ag thg previous sec tion, We 

restrict to 0 < si < \ and set a = (1 ~*^ , b = A. 

The typical term ~ j (1 — 5 0 ) k+j+1 produces the integral 

rh rb d5o 


dsi 


x 


Jo y/b — S 2 

S2 5 2a+k-l-j _ § jk+j +1 


y/a - hi 


dhx + 5 2 2 “ +fc - 1 - i 


( 5 i - S 2 ) 


k+j+l 


dsi 


S 2a+k-l-3 d ^ ,6 


2 

k+j+l 


'5 2 


y/a - hi 


-d<5i 


dh 2 


'0 


y/a — hi 


>h 


y/b — S 2 


rb +2cx+k—l — 


+ / dsi 
Jo 

= B 


j d5 2 [ a {8 1 ~5 2 ) 


k+j+l 


( k + j + 2, - 


1 




'$2 


y/a — 5\ 


-d5i 


dsi 


x 




k-\-j-\-3/2 


y/a — 5] 


+ / {a —8 2 ) 


k+j+ 3/2 o 2 


c2a:+fc—1 — 




yfb^-W 2 


(k + j + 1)! f 2 


© 


dsi 


fc+j+2 I '° 


X / V 


2a+k—l—j 


(b-v) t+1+I + 


dn 


(& +j + 1 )! 

© 


/ (2a + k — 1 — j, A; + j + 1). 


fc+J+2 


(Note 5 —.) The proof and statement of the I (in, n) formula are in 

V 2 / rn+1 

Prop lV.ll Then 


Pr {det / PT < 0 j = 2c ay kB (a, k + 1) 


a—1 


. (1 — a). (k + j + 1)! 

x ( -1 ) 7TTNT- m - 1 ( 2 ct + k-l-j,k + j + l) 

j =0 ^ )j+i ( 2 ) fc+y+2 


By formula (IV.ID 


I (2a + k — 1 - j, k + j + 1) = 


(2a + A; — 1 — j)! 2 (2 


a 


V 2/ k+j 


24* + « + 3 (4a + 3 k-j + l){(l) 
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We combine the various terms to arrive at: 


a—1 


2T (2a + 2 k + 2 ) 2 T (2a + k - jf 


Pr {det £ PT < 0} = - , 

1 J “J V (a + k + l ) 3 T (cc — j) T (4a + 3k — j + 2) 

T (2a + 2k + 3 ) 2 


2r (a + k + 2) T (4a + 3 k + 2) 

a—1 


X 


Et— 


-i) 3 


3=0 


(a — 1 )! (a + k + 1 ) 


(a + k + l) 2 _ x _ • (—1 — 3 k — 4a) ■; 


the j-sum is a polynomial of degree 2 a —3 — j in k (the factor (a + k + 1 ) obviously cancels 
when j < a— 2, and when j = a—1 the last factor equals —3 (a + k + 1) (—1 — 3k — 4a) Q _ 2 ). 
This phenomenon is qualitatively similar to our results for the general 4x4 situation. In 
particular for a = 1 


Pr {det £ pt < 0} 


2r (2k + 4 ) 2 
T (k + 2) r (3k + 6 ) 


2 " +7 (i)L 

Q3fe+5 f n ( 2 ) 

V 3/ k+2 V 3 / k+2 


and for a = 2 


Pr {det £ pt < 0} 


2 (k + 6 ) T ( 2 k + 6) 2 

3 T (k + 3) T (3k + 9)' 


A. Probability of det£ PT > det£ 


Next we compute Pr {det £ pp > det £} for a = 1 and the measure (det C,) k du a . As 
above, it is neater to work with the complement. From equation ([2]), det£ PT — det£ = 
(<!>i — <5 2 ) (s 5^2 — S 4 C>i). To make det £ PT — det £ < 0 either 5\ > 62 and S 4 > ^-55 or hi < f > 2 
and s 5 > |^s 4 . The first two integrations for > S 2 are 


(S 1 S 2 


\fc+i 


(l-s 5 ) k ds 5 (1 — s 4 ) fc ds 4 


r s 2 


s 5 


(SiS 2 ) 

fc+i 

k + 1 

rfc+i 

°2 

f 1 

k + 1 

Jo 

1 

k +1 

E 

k + 1 

3 =0 

1 

fc+i 



(1 - s 5 ) k ( 1 - J-S 5 


fc +1 


dsn 


(1 — S 5 ) 1 " (^1 — 82 + 82 (1 — Ss)) fc+1 ds 5 


k + 1 


k + 1 


(Si - s 2 y sf 2 


j <-2k+2-j 


(! - s 5 ) 


2 k+ 1 ~ j ds, 


3=0 


k + l\ (5 1 -8 2 ) j 5 2 2 k+2 - j 
2 k + 2 — j 
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similarly for <h < 82 the value of the integral is Y^ k j =o ( fc + 1 ) P 2 -• Proceeding as 

before and with the same notations carry out the remaining integrations on the term with 
fixed j, for 0 < j < k + 1. Assume 0 < si < then we obtain 


dsi 


d<5 2 


b5 2 r2fc+2 —j 




2 J 0 


St ^ 3 (h - Si ) 3 

y/a — hi 


dh! + Sf K ^~ J 


2fc+2 -i r & - w 


' $2 


\/a — 81 


■d<5i 


„ b 6 2 k+ 2 - j d6i [b fa-Si ) 3 <162 f b S'f/ + ' z - J dS 2 f a ( 8 ,- 82 ) 


b c2k+2—j 


dsi 

Jo 

= B [j + 1 , 


+ 


y/ a — 8 1 Js 1 \/b — S 2 

1 \ rb ( r b 

' As A (b-S^ 6 -! 


'0 y/b — S 2 Js 2 V a — hi 


-dhi 


2 k+2-j^c 


y/a — 8 , J 0 


+ / (a - 82 ) 


j+ 1/2 °2 


c2fc+2- 


J dh 2 


y/b ~ 82 


V- 


(h 


dsi / 


2 k+2—j 


3 +1 




dn 


= -^— J(2A: + 2 - j,j) 

UJj+i 


Thus 


fc+i 


/& + 1\ 


Pr{det « PT ^ det d = 2 v*E(d j (t + 1)(2fc + 2 _.,- )( . )j+ /< 2fc + 2 -m 


_ n2k-l (Dfc (Dfc ^ ®) 

j= [ 


3 


^+ 1 ) ! (l) 2fc+ l^(- 2 ^ 2 ),(- 2 ^ 1 )/ 


after simplication. The sum is a truncated 3 -F 2 series. 

With more techical details, we can determine Pr {det / PT < det /} for the measure 
(det /) k dv a for ct, k = 1, 2, 3,.... Proceeding as for a = 1 we start with the integral 


(<M 2 ) 


a+k 


S 5 (1 - s 5 ) ds 5 / S 4 (1 - s 4 ) ds 4 , 

J ss&o 


where h 0 = ^ and 0 < 8 2 < hi; and there is a similar expression when 8 , < 8 2 - Use formula 
(J5]) (with <5 0 replaced by s 5 h 0 ) to get the value 


a— 1 


4 ' 


‘"-’■S'-"'Afe 


82 1 ^ s 5 1 J (hi — h 2 + (1 — S5) h 2 ) fc+,7 + 1 dss 


(1 a)r Yt{ k + l+j )B{2 a -i-l,2k + j-i + 2)(S 1 -6 2 ) i ^ 


_ _ _ i k+j+l 

§ ( d) (fc+ip+i u v * 

k+Ot. 

^2 {Si -8 2 ) i 8 2 2 k+2a ~ i a(k,a,i), 


i =0 


14 
































where 


a—1 


cr(k,a,i) = ^ 


i 7 =max( 0 ,z— k— 1) 


k\ (a — 1)! (2a — 2 — j)\ (2 k + j — i + 1)! 
z! (a — 1 — j)! (k + 1 + j — i)\ (2k + 2 a — i)\ 


For 0 < i < k 


k\ (2k + 1 — z)! (2a — 2)! (1 - a). (2k + 2 - i) j 

(7 Q^, 1) .j , 1 ; \| y 07 , <->^ i \! / v 


z! (fc + 1 - z)! (2k + 2a - 1)! ^ (2 - 2a), (fc + 2 - i) ’ 


and for A; + l<z<& + a 


cr (fc, a, z) = < a 


j=o 

k + a x ^ _2 
a 

fc+a /y X \i c2/c+2a —i 



k + a 
z 


Similarly, the integral equals (<5 2 — “ V(fc,a,z) when < 5 2 . With the 

same steps as above, we obtain 

k+a .| 

Pr {det £ PT < det £} = 2 c a ^ ° (^><*) *) tjt —I (2fc + 2a — z, z). 

i=0 \2/j 

Some simplification may be possible. 


'»+i 


V. INTEGRAL FORMULAS 


In this section we compute closed expressions for 


I( m ,n):= I f dsj^ v"‘\(a-vY fj^+ib-vY^^Av, 

for m, n = 0,1, 2, 3,.. . and for (m, rz) = (2a, 0) with a > 0; where a = (^) 2 , & = 
that 0 < b < a < The auxiliary formula 

S n ) (~ m h ( n + 1 )» ^ ~ j 1 

' ’ 2-^ n\ (m I n I 9 v2j 


2> SO 


^z!(m + n + 2),.^ j! (|) * + .? + § 


= 2 


2m+2n 


m! (m + n)! (m + n 


is proved in 


20 


n \ ( n + 2m+ l)!Q) m+n+1 


, Prop. 2], 


Proposition V.l For m, n = 0,1,2, 3,... 


I (m, n ) — j ds J v m ^ (a — v) n ^ ^—- + (6 — v) r 


a — v 


b — v 


a — v 


dn 


= 2 


—4m—4n—5 


= 2 


—2m—2n—3 


B (m + 1, n + 2) S ( m , n + 1) 
(m!) 2 (m + n + 1)! (f) n+1 


( 2 m + n + 2 )! (|) m+n+2 
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Proof Set z = 1 — 2s, u 2 = lj ~z L , a = p (1 + z) 2 , b = -b (1 — z) 2 then 


b — v = 


u 2 z 


4(1 -v?)' 

2 „,2 


a — v = 


4 (1 — -u 2 ) ’ 

(1 - z) z - u 2 (1 + z) 2 
v =-——-rr-, du = 


uz 


16(1 -u 2 ) 


2(1 -u 2 ) 


rd u, 


0 < 2 < 1,0 < u < 


1 -z 
1 + z' 


For changing the order of integration note that 0 < z < 1 and 0 < u < is equivalent to 
0 < u < 1 and 0 < z<\^. Thus 

— — — — 1+U 


(a - vf J+ {b- v) n \ / h — \ dsdu 


b — v 


u 2 z 


u 2 + 


a — v 

z 


zdzdu 


_ _ 2 ~ 2 n- 2 z n+l 


4(1 -u 2 )J \4{l-u 2 )J J 4(1 — u 2 Y 

u 2n+2 + 1 


(1-u 2 ) 


2 \n+2 


dzdu. 


Then introduce the new variable y: 


1 — u z n+1 . y" . 

z = T—y > t: -= - - . 2n i 4 dy, 


,n +1 


l + u*’ (1-U 2 ) n+2 (1 + u) 

v = 15 (1 - y) (* - (rHO y 

0 < u < 1,0 < y < 1, 



and the integral is transformed to 


*1 „.2n+2 


)—4m—2n—2 


u zn+z + 1 


d u I y n+1 (1 - y) m 1 - 


/o (1 + u) 

_ 2 ~4m— 2n— 2 


2n+4 


B (m + 1, n + 2) ^ — 


(“Hi + 2 ); 


i=0 



( 6 ) 


2i 


du, 


by use of the binomial expansion and 
/■i 

,n+l 


(" + 2 )i 


I y (1 - y) m y l &y = B (m + l,n + 2 + i) = B (m + l,n + 2) ^ + ^ + ^ 
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T 1 — u 1 — w , , 

Let w =-, u =-, then cm = 

1 + u 1 + w 


2 1 

idw, —— = \ (1 + w), and 


(l + tu)“ 1 + u 


n M 2n+ 2 + 1 n_ u 
, (1 + u) 2n+A V 1 + u 


2 i 


du = 2 _2n_3 


{(1 - w) 2n+2 + (1 + w) 2n+2 } w 2i dw 


„1 n+l 


= 2 


—2 n —3 


2 ^ f 2 ’\ + 2 )w 2 J +2 ’dw 


'« V 2 i > 


= 0 — 2 ) 1 — 3 v { — ~ 1)j + Mi 

- + i+)++i+) - 


'3 


Thus the integral equals 


n-4m-4n-5 n / , i , o') V"' ( _m )i ( n + 2 )i V'' ^ 71 ( U + 

2 B ( m+ l,„ + 2)g. |(m + n + 3) X Jl( . ) (i + j + i) 


= 2 


—4m—4n—5 


.B (m + 1, n + 2) S'(m, n + 1). 


The proof is completed by simplifying F? (m + 1, n + 2) S (m, n + l). 


Proposition V.2 For a > 0, 


I (2a, 0) 


7 T T(2a+1) 2 

2 6+8 " T (2a + |) T (2a + §)' 


Proof Follow the previous proof from the beginning to equation ((HD which is replaced by 



= 2~ 8a ~ 2 B (2,2a+1)^2 

i =0 


(-2a).(2)i 

i! (2a + 3) j 


C )< 2 - 1 

o (1+u) 4 


1 — u 

1 + M 


2 i 


du. 


The binomial series terminates when 2a e Z + ; otherwise from the asymptotic property of 
the gamma function it follows that (with some fixed m € Z + and m — 2a > 0) 


(-2a) 


m+n 


(m + n)! 


(~2a) m 

(m + n)! 
(-2a) 


?n 


2 “)„ = 


m!T (m — 2a) 


(—2a) m T(m-2a + n) 

m!T (m — 2a) T (m + 1 + n) 

n~ 2a ~1 , (n —>■ oo) , 
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so by the comparison test the series converges for 0 < y < 1 and 0 < u < 1 provided 
a > 0, and can be integrated term-by-term. Also B (2, 2a + 1) = (0a+1) 1 ( . ?a+i - l) With the same 
argument as in the previous proof, with n = 0 we obtain 


u 2 + 1 


(1 + uf \l + u 


u 


2 i 


d u = 


i + 


+ 


* + | 


0 + 1 ) 


4 () - i) (l - |) 


a + 1)! (l). (2), (h, 
«(i) 


i+2 


Then 


/ ( 2 a, 0 ) = 


)—8q—2 


• ^ o ! 


3(l)i(!)f 

- 2 a) i ( 2 ) i ( 2 ) j (l) < 


3 (2a+1) (2a+ 2 ) ^ i! (2a + 3),(l), (f)/ 

This is^a very-well-poised 4 T 3 series which is summable using the Rogers-Dougall formula 
(see 


lg 16.4.9]): 

5F4 


a, | + 1 , b, c, d 


; 1 


, a — 6 + 1,a — c + 1, a — d + V 
T(a — 6+1) T (a — c+l)T (a — d+l)r (a — 6 — c — d+1) 


T(a + l)T(a — b — c+l)T(a — b — d+l)T(a — c — d + 1) ’ 

the sum converges if one of b , c,d is a negative integer and it terminates or b + c + d — a < 1. 
The formula applies to our sum with a = 2, b = —2a, c—\,d — | (that is, d = a — d + 1), 
and b + c + d — a = —2a < 1. The result is 

2- 8q -2 r (2a + 3) T (|) T (|) T (2a + 1) 


I (2a, 0 ) = 


3 (2a + 1) (2a + 2) T (3) T (2a + §) T (2a + §) T (1) 
T(2a + 3)T(2a + 1) 


2—8o—6 


7 r 


(2a + 1 ) (2a + 2 ) T (2a + §) T ( 2 a + §) 
T( 2 a + l ) 2 * * 


7 r 


2 6 + 8 « T ( 2 a + |) T ( 2 a + §) 

Of course the two formulas agree when 2a G Z> 0 . 
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